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Abstract 

O' 

Qv^ ■ Under the framework of the Bethe-Salpeter (B.S.) wave functions and the 

i-Q \ Mandelstam formaUsm as well, to make "instantaneous approximation" to 

f^ • a transition matrix element (a current operator sandwiched between two 

(— I \ bound-states of double heavy quarks) is described. By taking the typical 



concerned decays as examples, such as Be meson decaying to J/V' + i}^)-, the 
advantages of the approach and its limitations are illustrated. Finally, po- 
tential applications to various processes for possible double heavy flavoured 
systems, such as those of (Q'Q) and (Q'Q) {mQ^niQi » ^qcd), are dis- 
cussed. 
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The Standard Model (S.M.), as well-known, has made great successful. The 
flavour decays i.e. the weak decays, of the fundamental fermions (quarks and lep- 
tons) are supposed to be understood very well in the framework of the S.M. i.e. the 
vector boson exchange and the CKM (Cabibbo-Kobayashi-Maskawa) mixing ma- 
trix structure. However, except the leptons and the undiscovered top quark, all the 
known quarks always decay as constituenets of hadrons, once they are produced. 
This is due to the fact that colour has confinement nature and the quarks carry 
colour explicitly, thus a quark always 'fragment' into hadrons before decaying once 
it is produced. (Only the undiscovered top quark may be an except, as its mass is 
very great hence its lifetime is so short that there is no time to fragment into hadrons 
before decaying.) Thus we always see flavour decays through hadrons in experiments 
and we have to deal with the effects of the confinement when studying the flavour 
decays of quarks. However, the confinement effects are of a nonperturbative nature 
so very difficult to compute them reliablely. Hence the weak decay problem has a 
very long story, which had started before the S.M. was established. Nevertheless the 
progresses on the problem have always being achieved steadly, especially the recent 
achievements of the Heavy Effective Theory (HET) for one (or two) light and one 
heavy flavour hadrons: Qq or Qqq {rriQ >> rriq) and their antiparticles. Due to the 
authors of HET'^^^1, a great progress has been achieved in understanding the weak 
decay of heavy flavours under the approximation rriQ -^ oo and with the corrections 
of the order l/mq and higher. The authors note at very beginning when they pro- 
posed the approach that one of the applicable conditions of their formalism is that 
the mass of the heavy quark in the concerned hadrons should be greater than Aqcd 
and much greater than that of the light one. 

In the meantime as emphasized by many authors[7-ll], the B^ meson will play 
a special role for understanding the heavy flavour decays and some progresses have 
achieved too. 



The Be meson is the ground state of the bound states for the system of the 
heavy quark-antiquark pair c and b. Both of the quarks b and c are heavy quarks 
i.e. rriQ <^ Aqcd, Q = b,c. Though both may be considered being nonrelativistic, 
their masses are compatible, thus the formahsm^^"^] may not be very appropriate, at 
least careful examination are needed and considerable modifications are expected. 

Of the three generations in the S.M., besides the well-studied systems of the bb 
and the cc, the quark-antiquark c and b and their antiparticles are the only ones 
which still have a long enough lifetime to form a bound state before one of them 
decaying. In contrary, the rest possible heavy quark pair systems must contain a 
top or anti-top quark at least. They will not have chance to form a bound state. It 
is due to a great mass of the top quark ( the top mass rrit greater than lAOGeV^^'^\ 
though the top has not discovered yet ) that causes a great phase space so a very 
short lifetime even in the case only weak decays are allowed, hence there is no time 
to form a bound state before the top decaying. The c and b system so interesting 
is attributed to the fact that it is the only possible one which may make mesons 
and the mesons have not been discovered yet. Besides it plays a special role in 
understanding the heavy fiavour decays. The Be meson, the ground state of the c 
and b system, not as that of the c(or b) and c{b) systems, carries flavors explicitly 
so no strong decay at all i.e. the ground state is stable to the strong interaction. It 
has a great advantage for studying the weak decay and the bound state effects of a 
heavy quark-antiquark system, even for the later only, still better than those of the 
J/ip and T families. Indeed, the ground state, the Be meson, has many sizable and 
various weak decay channels, so that it becomes possible to compare the results of 
experiments and theories on the weak interaction effects and the binding effects of 
a heavy quark pair system widely. 

Usually the above two kinds of effects in the meson decays are divided into 
two categories from another aspect for convenience: the short distance effects i.e. 



the weak interaction responding for the flavour change and its perturbative QCD 
corrections; and the long distance ones i.e. the binding effects for the components. 
In general speaking, although they are entangled, the second ones, being of non- 
perturbative nature, are very difficult to deal with aliably, while the former ones, 
being perturbative nature, are relatively easy to treat by means of the Standard 
Model of electro-weak interaction and perturbative QCD up to a desired accuracy. 
However, for the non-perturbative binding effects of a heavy quark system, great 
progresses have been achieved under the framework of the QCD inspired potential 
model, but in some phenomenological sense, especially for the J /ip and T families. 
Thus we may consult the potential model to treat the binding effects for the c and 
h system too. According to the QCD inspired potential model, the Be meson is 
very similar to the two families, of J/ip and T, but has a different reduce mass 
only. Thus the Shrodinger equation with the same QCD inspired potential and 
a corresponding different reduce mass, should be worked very well to depict the 
bending effects of the c and h system. The Shrodinger equation, as known being 
in a sense of nonrelativistics, hence to the static properties of the system, is good 
enough, while to the decays, especially those with a great recoil momentum may be 
problematic when using the solution(s) as the wave function(s) directly. To find a 
way to use the Schodinger wave function(s) properly into the decays of the double 
heavy quark systems is the aim of the paper. 

To our knowledge, the Mandelstam formalism[13] with the Bethe-Salpeter (B.S.) 
wave functions is one of the best approaches to compute a matrix element sand- 
wiched between two single bound states, however which is very often to be attributed 
at last when calculating a whole decay matrix element. Due to Salpeter[14], the re- 
lation of a B.S. wave function to the Shrodinger one was established quite long time 
ago by the so-called instantaneous approximation. Therefore it seems no problem 
for the present purpose that the Mandelstam formalism might be used with the help 



of the Shrodinger wave functions by solving the reduced and corresponding nonrela- 
tivistic equation. Nevertheless, in the Salpeter's approach, the approximation taken 
at the centre mass system of the bound state is implicated. If only to apply the 
Shrodinger wave functions to calculating those decay matrix elements sandwiched 
between two states, but of which only one bound state is involved, it is straightfor- 
ward relatively, provided one calculates such a matrix element at the centre mass 
system of the bound state. However now for the most important decays of the Be 
meson, the problem is not so straightforward because the decay matrix element is 
attributed to calculate a weak current matrix element with quite a great momentum 
recoil, the current operator sandwiched between two bound state. In this paper we 
are generalizing the instantaneous approximation to suit the interesting cases. 

To be the start point, we suppose that the transition matrix elements of the Be 
meson, sandwiched between two single bound states, ( the form factors or the decay 
current matrix elements ) may be computed in a desired accuracy with the help of 
the Mandelstam formalism no matter how great the momentum recoil, provided that 
the one-particle irreducible green function, responsible for the matrix element in the 
formalism, matches to the trancated kernel of the B.S. equation well, as required 
by the formalism'^^^ . It is because the formalism is fully relativistic. Now the prob- 
lem is attributed to having reliable B.S. wave functions for the quark bound states, 
appearing in the formalism. In general as for hadrons the key problem to obtain 
reliable B.S. wave functions is to establish the B.S. equation with a reliable kernel 
and to find a suitable method to solve the equation. However the nonperturbative 
nature of the binding effects between the quarks from the QCD confinement makes 
the problem very difficult. For those matrix elements sandwiched between two dou- 
ble heavy quark bound state(s) such as the B^ meson and J/ip ( or T ) etc. the 
B.S. wave functions occurring in the formalism can be acquired reliablly by means 
of the instantaneous approximation, i.e. the B.S. kernel for the double heavy quark 



systems may be related to the QCD inspired potential of the Schodinger equation 
esteblished in a potential model framework, so the B.S.'s and Schrodinger's wave 
functions are related too through the approximation. Whereas one could not set the 
initial and final states in a reference frame both at rest for a real decay process, espe- 
cially when a matrix element with quite a great momentum recoil is considered, thus 
the relation between the wave functions of the potential model and the B.S. ones, 
established only at a rest frame (C.M.S.) as indicated by the original instantaneous 
approximation, still cannot be applied straightforward. Here we constrain ourselves 
to consider the transitions between the states of double heavy quark systems. The 
strategy, we adopt is to generalize the instantaneous approximation to apply to the 
whole matrix elements no matter how great the momentum recoil. We note here that 
as the wave functions are computed out only numerically in the potential model, one 
cannot use the method by Lorentz boosting it directly to pursue the present purpose 
reliably. Our approach here, essentially to say, is to calculate the matrix elements 
under the Mandelstam formalism and to apply the 'instantaneous approximation' 
to the matrix elements in a whole, instead of that to the B.S. equation only, as done 
by Salpeter'^^1. One will see that with our approach the calculations on the matrix 
elements under the Mandelstam formalism and the B.S. wave functions extracted 
from potential models as well, are reliable even in the cases of the transitions with 
quite a great momentum recoil. As previously discussed, we develop the approach 
with the motivation from the Be meson physics, thus we will take the decays such 
as the Be meson to J/ip + (lu) as examples to illustrate it in this paper. 

To calculate the exclusive weak decays of the Be meson, one needs to evaluate 
the hadronic matrix elements, i.e. the weak current operator sandwiched between 
the initial state of the B^ meson and the concerned hadronic finial state. We restrain 
ourselves to evaluate them in the simplest cases, i.e. only those decays in their final 
state having only one hadron for semileptonic ones, but two hadrons for nonleptonic 



ones. In these cases, one may attribute the problem to evaluating a matrix element 
of the weak current operator sandwiched by two of single-hadron states ( for non- 
leptonic decays it is due to the factorization assumption of calculating the decay 
amplitude ). 

With the notation of a weak charged current J^ = V^ — A^, where V^, A^ are 
the vector and the axial vector current respectively, the matrix elements are related 
to the form factors'^^'^^^ as 

< Pip')\V,\B,ip) >= Uip + p'), + f-{p- p')„ 

< V{p',e*)\V^\B,{p) >= ige^,p„e*\p + p')P{p-p'Y, (1) 

< V{p', e*)\A^\B,{p) >= /e* + a+(e* ■ p){p + p')^ + a_(e* ■ p){p - p')^. 

where p, p' are the momenta of the B,. and the outgoing hadron respectively, P and 
V denote the pseudoscalar and the vector mesons respectively, e is the polarization 
vector of the vector meson. The form factors are functions of Lorentz invariant 
variable r^ = [p — p'^. 

In the literature one may find two kinds of approaches to calculate these form 
factors. One of them is BSW model'^^', in which the authors calculated the form 
factors at the maximum recoil r^ = by means of the wave functions defined at the 
light cone system under the quark model framework, and then to extrapolate the 
result to all values of r^ by assuming the form factors dominated by a proper pole of 
the nearest ones. The other is IGSW modelt^^'. The authors of Ref. [17] calculated 
the form factors by using the wave functions of the quark-model ( "mock meson" ) 
which treats the hadrons as a nonrelativistic object. As argued by the authors, the 
approach is exactly valid in the limit of weak binding and at the point of zero recoil. 
However, in the cases with a large recoil, it is problematic. For instance, for the 
decay Be -^ J/ip + p, i.e. the example we are considering, although the initial state 
Be and the final state J/ip both are of weak binding, the recoil of the decay is not 
small so the IGSW model may not work very well. 
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However, inside the Be meson, both the b and c are heavy so may be con- 
sidered being nonrelativistic objects, but their masses are comparable, thus the 
formahsm^^^®] may not be very appropriate, at least careful examination and con- 
siderable modification are expected. Furthermore, even though their approach is 
not suitable for the B^ meson, to establish a link between the universal Isgur-Wise 
function ^{v ■ v') and the nonrelativistic wave function of the heavy meson obtained 
by the potential model would be still very interesting, where v and v' are the four- 
velocity vectors of the meson in the initial and final states respectively. 

Now let us proceed to write down the matrix element with the help of the Man- 
delstam formalism, so to describe the approach explicitly. It is known that the B-S 
equation of a fermion-antifermion bound state takes the following form: 

/■ d^k 

{j)i - mi) Xp{q) ih + m2) = i J —^ V{p, k, q) Xp{k) (2) 

where pi, and p2 are the momenta of the constituent particles 1 and 2 respectively. 

They can be expressed in terms of the total and the relative momenta p and q as 

mi 
Pi = aip + q, ai 



m-i + mo 

(3) 
m2 

P2 = Oi2p-q, a2 = ■ , 

m,i + m,2 

V{p, k, q) is the interaction kernel. It is well known that the B-S wave function Xp{q) 
satisfies the normalization condition: 

-§^^^MU<l)^^[Si\pi)S2\p2)6\q-q'^^^ = 2tpo. (4) 

As for the decays caused by the quark gi of the meson, according to the Mandel- 
stam formalism[13] and the spectator mechanism shown in Fig.l, the weak current 
matrix element involving one hadron in the initial state and one in the final state 
respectively, may be expressed in terms of the B-S wave functions: 

nr)=^J-0^Mx'pi^')^'tXp{q){h + rn2)], (5) 
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where Xpio)) Xp'W) ^-^e the B-S wave functions of the initial state and the final state 
with the total momenta p, p' and the relative momenta q, q' respectively; pi, mi, p[, 
m[,and p2, m2 are the momenta and the masses of the decay quark, the final one 
and the spectator respectively; T^ is the weak interaction vertex and to the lowest 
order, Ti has the form of 7^(1 — 75) for the charged current. 

As pointed out above, the B-S wave function Xpil) of the heavy quark pair 
system under the instantaneous approximation, can be evaluated by solving the 
corresponding Schrodinger equation with a QCD inspired potential. Usually, to 
make the nonrelativistic instantaneous approximation for the B-S equation is implied 
in the rest frame of the concerned bound state. 

Now let us outline the the approximation but in a form, which may be used 
further. 

The instantaneous approximation is to carry through the integration over the go 
component for the B-S equation Eq.(2) firstly, when the kernel at the rest frame has 
a simple form 

l^(p,A;,g)~V(|k-q|), (6) 

then as a result, the B-S equation Eq.(2) is deduced into a three-dimensional equa- 
tion, i.e., the Schrodinger equation in a momentum representationt^^l . 

However, as pointed out above, when the decays with sizable recoil are concerned, 
it still is not sufficient to have the relation between the wave functions in the rest 
frame only, because in the Mandelstam formalism Eq.(5), the start point of our 
approach, at least one of the B-S wave functions in a Lorentz covariant form in a 
moving frame is required. To pursue the purpose, we take the way to make the 
instantaneous approximation to the whole matrix element itself instead of only the 
wave functions. 

As the first step, we need to divide the relative momentum q of the B-S equation 
into two covariant parts firstly: qp\\ and qp±, a parallel one and an orthogonal one 



to the centre mass momentum p respectively, i.e., 

q' = <ii + q'pL (7) 

where q^u = -jj^p'^', qp± = Q^ — C[p\\- Correspondingly, we have two Lorentz invariant 

p 

variables: 

_ P^ 

Mp (8) 

(IpT = ^ql - q^ = ^-qp±- 
In the rest frame of the meson, i.e. p = 0, they turn back to the usual component 
go and |q| respectively. 

Now in terms of these variables, the covariant form of the wave function can be 
obtained. The volume element of the relative momentum k can be written in an 
invariant form: 

d'^k = dkpkp-ydkpjdsdcl), (9) 

where is the azimuthal angle, s = -^^ . The interaction kernel can be 

kpjqpT 

denoted as: 

V{\k-q\) = Vikp^,s,qp^), (10) 

which is independent of kp and Qp. Let us introduce the notation for convenience: 

dqp^ 



'Ppiqp±) = ^j^Xp{qp\\,qp±), 



v{qp±) = J — (27r)3 '^(V.g.gp^)yp(^p±) 



;ii) 



thus the B-S equation can be rewritten in short as: 

Xp{qp\\,qp±) = Si{pi)T]{qpA_)S2{p2), (12) 

where 5*1(^1) and 6*2(^2) are the propagators of the constituent particles and they 
can be decomposed as: 

gf^\ ^ Kpi.qp±) ^ Kp{.qpA-) .^gx 

J{i)qp + ctjAf - uoip + it J{i)qp + onM + uoip - ie ' 



with 



UJ. 



ip 



mf + ql-j 



Kil: 



1, 2 and J(i) 



ip\.'dp±) 



2u. 



ip 






(14) 



Here A^{qp±) satisfies the following relations: 



Kp{(lp±) + Kp{Qp± 



M' 



^ip{<lp±)j^^ip{<lp±) - Kp{ 



(15) 



A-tp{qp±)YJ^fpi<lp±) = 0. 



M 



Thus, Ajp(gp^) may be referred to as p-projection operators ( p is the momentum of 
the bound state ) while in the rest frame they correspond to the energy projection 
operators. 



If defining (Pp^{qp±) as 



^p^iqp±) = A-ipiqp±){TMlp±)^fpi(lp±) 



M'' 



M' 



(16) 



where the upper index c denotes the charge conjugation, under the notation: 



A5(gp±) = A|(gp±), 



and integrating over qp on both sides of Eq.(12), we obtain 

{M -uJip-uJ2p)'^p^{qpA_) = Aip{qp±)Vp{qp±)^2p{qp±), 
{M + uJip + uj2p)(Pp'{qp±) = Aip{qp±)Vp{qp±)Kp{qp±), 
Vp~{qp±) =Vp^{qp±) = 0. 

The normalization condition becomes into the following covariant form: 



f qljdqpT 
J 27r2 



<+(gp±)^<+(?p±)^-v5p {qpi-)j^Vp iqp±)j^ 



(17) 



2M. (1^ 
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Now for the usage later on, let us introduce two 3-momenta pi and p2, 

(19) 



Pi = ^p + qp±> 



^2p 

P2 = ^P - qpx. 

In the case of the weak binding as concerned here, the wave functions of the double 
heavy quark systems can be constructed approximately as follows: 



1 

Vp~~iqp±) = J2 /4^ ^ ■■VsiPl)Us'ip2)(l);iqpT)xi': 



(20) 



where UsiPi), Vs'{pi) (i=l,2) are the Dirac spinors of free particles with masses rrii 
and momenta p,; Xss' i^ ^^e Clebsch-Gordan coefficients that make s' and s couple 
to A; and 0^(q'pT) is the scalar part of the wave function. 

In the present case, owing to weak binding, the ip^ (q'p±) is a small component 
and can be ignored. In fact, if the kernel is of scalar and/or vector, the ip^ {qp±) 
is in the order of (v/c)^ to 'Pp^'^ {qp±y^^^ ■ Furthermore, if ignoring the components 
proportional to the qp± in the spinor structure due to the nonrelativistic nature, 
(/9p++(gp_L) can be simplified: 

^i^Hqp±) = ^^(«75 + /3/)0(gpT) (21) 

where a = 1, j3 = for an ^5*0 state and a = 0, /3 = 1 for a ^5*1 state, while the 
"radius" wave function (piqpj) satisfies the following Schrodinger equation: 

|^0(gpT) + J ""^^ 'I Vis, V, qp±)<PikpT) = E<P{qpj), (22) 

with the reduced mass fi = of the system. 

rrii + m2 

Thus, we have established the relation in covariant form between the B-S wave 
function and the solution of the Schrodinger equation with the instantaneous ap- 
proximation. 
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Now we start the procedure to make the instantaneous approximation to the 
whole matrix element itself. The weak current matrix element, by integrating over 
the go component of Eq.(lO) with the method of contour integration. 

As asserted above, the negative energy parts of the wave functions are very 
small in weak binding cases so that we can ignore their contributions in the lowest 
order approximation. In the present case, we are concerning weak-binding bound 
states only, thus those contributions from the negative energy parts will be ignored. 
Putting Eqs.(12), (13) into Eq.(5), we have the matrix element: 



1 (r) - I ^ \n' (a' ) ^^A%'^^ y^ 

q , + aiM - urp, + le 



(27r)4 ['P^" 

Kp{(lp±) A^p(-gpx) 



Qp + aiM + uip + ie Qp - a2M + uj2p - ie 



VpiQ', 



'p±) 



In the bracket of the integrand, there are three poles at points Oj in the complex- 
go plane: 



ai = —aiM + cji — ie, 
02 = a2M - UJ2 + ^e, 



(24) 



a[ = a2M - Ey (r + q)2 + m[^ - ie, 
and two branch cuts starting at the branch points: 

go~m2±i^-, (25) 

7 

Irl 
with 7 = — and all the terms relevant to q are ignored, due to the fact u'l , = 

^Jq'p^± + rn^- 

Now the problem becomes how to carry out the integration over the go component 
on the right side of Eq.(23) or, in another words, how to treat the cuts of the 
integrand if contour integration method is adopted. Wereas here we will treat the 
branch cuts in the integrand approximately, by expanding uj'ip, as follows: 

'2 



UJ 



J 



(1. 



qph + m? = m[ + ^ + ... . (26) 
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In a weak binding limit it is quite good approximation. According to Cauchy's 
theorem, the integration of a closing contour on the upper half plane of the complex- 
go for the matrix element lfj,{r), is just that summing up all poles' residues. However, 
as the pole 02 on the upper half plane is very close to the pole Oi on the lower half 
plane, i.e. the distance 

A = ai - aa ~ M - ml - m2 + -^ (27) 

2/ii 

is small, the value of the integration is dominated by the residue of the pole 02 only. 
The contribution from the pole a' is not important i.e. may be ignored at all at the 
concerned accurate level.0. Therefore, we obtain: 

J {2Tcy '' ^ {qp, + a[M'-uj[p,){M-Uip-uj2p) 
It is easy to prove the following relations 



^L 



(29) 



A^^(-gpx) = ^A^2?(-9;'±)7oA2?(-gpJ , 

and with Eq.(17) from Eq.(28) we obtain the required equation: 

«^) - / nip^" I "J^f'-'r-^r (*J^ 1 J. (30) 

where qp,j_, u[p,, uj'2p/, ujip, uj2p are as expressed in Eqs.(8), (14) and (26) i.e. 

/ / '2 \ ^ 



^iv = ^^ + ^h 



, _ E'uj2p + r ■ q 

^2p' - J^, , 



/ I '2 2 

Vt — \l^2p' ~~ ^2- 



(31) 



^ In fact, some poles ( even cuts ) may be induced into the matrix element through fj'i{<l',j_), 
however as the similar reason as here, they would not contribute so substantially as that of the 
pole 02 to the final results in all the cases of weak binding. 

13 



We should note here that based on the adopted extra approximations such as 

Irl 
Eq.(25), the Eq.(30) is vahd only with not too large recoils, i.e. 7 = — — < 1, 

however most of the decays of double heavy quark systems such as our concerning 

processes Be -^ J/ip + X, Be ^ Bs + X etc. satisfy the condition well.0 Using the 

wave functions in the form of Eq.(20) for both of the initial state and the final state, 

it follows that 



^f^ir) = / ^ [ui(p'i)r^M„(pi) ] ^p{q'p>T)'^p{qpT)Xis'Xs'r 



[zn) ^uj2pUJpUJ^p, 

(32) 



where 



^ip' = JQp't + ^i 



Pi = (^i,q) 



(33) 



^Ip + ^2p I ^Ip + ^2p , 



and the normalization condition for the 'spectator' ( the antifermion with momentum 
P2 in Fig. 1 ), 

Vs'{P2)-r7Vv{p2) = 2U26s'V, (34) 

has been used. 

In order to compare with those results of Refs.[16,17], we compute out the form 
factors and it is very interesting that we may extract a similar 'universal' Isgur-Wise 
function at last. 

There is some arbitrariness in choosing the directions of the spins of the quarks 
when calculating the form factors. Thus we may use it to simplify the calculations. 
It is convenient to take direction of the spins orthogonal to the p and p' because the 
spins in this direction remains unchanged if a Lorentz boost along the p' direction 
is taken. 



^Here in order to show the approach emphasized in this paper briefly so as to see the inspire 
of the approach clearly we take the extra approximation, as matter of fact it may be weaken or 
improved in the practical usage. 
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In the Appendix of Ref. [7], a covariant formalism to calculate the creation of a 
pair of fermion-antifermion have been derived in the spirit of helicity amplitude^^^'. 
We will employ the method for the present problem. 

A similar formalism can be obtained for a fermion scattered by a virtue W^ 
boson i.e. the amplitudes with possible spin directions read as 



1,2 +2 

M,^ = L_-tr 
3,4 2 






(35) 



where 



L± = [-(pi-Pi±mim'J] 2, (36) 



and ki is an auxiliary and space-like vector {kf = — 1), which is orthogonal to the 
initial momentum pi and the final one p[. It is easy to see that the first equation of 
Eq.(25) describes the spin nonflip amplitude while the second describes a flip one. 
Both fermions are fully polarized along the ±ki directions. It should be noted here 
that in these formulae the relative phases of the spinors among those states with 
different polarizations have been fixed. 

To calculate the form factors from a current matrix element, we need to construct 
the spin wave functions of the individual quarks into definite spin ones to describe 
the initial and the final states of the matrix element respectively. In general, for an 
"^S'o state 

Xss' = ^(Ti - it), (37) 

and for an ^5*1 state, the spin structure corresponding to the three possible indepen- 
dent polarizations are 

x^sl' = ^(Ti + it), 

X'sl' = ^(TT + ii), (38) 

xS' = ^(TT - ii), 
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where three vectors ki, k2, and k^ orthogonal to each other, are used to denote the 
polarization directions of the ^5*1 state. 

Thus, for a transition P ^ P' -\- X , the amplitude reads 

M^ = L+itr[ if, + m;)7575(A + m^)T'^ ], (39) 

and, for a transition P ^' V + X, the corresponding amplitudes read 

Mf = L+itr[ (/i + m[)-fM^i + rm)r^ ], 

M^ = L_itr[ (/i + m[)M^, + rm)T^ ], (40) 

M^ = L_itr[ (/i + m;)75(|^i + m,)T>^ ], 

where Mi, M2 and M3 correspond to those of various projections of the polarizations 
of the final states. In fact, if letting 

k!l = h+L^e'^-'Pylk,,p^, 

(41) 

ks = -^m^Kp' ■ P'i)P'^ - (P' ■ P^Pi"^ 

correspond to the polarizations of M2 and M3, the amplitudes of the P ^ V + X 
transition can be written down in a compact form: 

M^ = L+^tr[ (/i + m[)-f,/{^, + rm)r^ ], (42) 



where 



(P ■ P') ■ 



e — e / /N Pp'_L, 



P^ '2 



C = ^+ - 1, (43) 

{p' ■ piW ■ p'l) 1 



M'2 L' 

(p-p') , 



2 



PP'±=P-^^^7^P. 
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After a straightforward calculation, the form factors may be formulated: 



i 



OJ, 



UJr, 



I ' 



/ 



a± 



MM'm'i 



' {P-P'i) 



21712 , ^ f iu[+iU2 jp-p') ^^ 



(44) 



M2m'i 



where 



C(l + 



MM'm[ 
jp-p'i) - 



M'2 



pz 



(P-P' 



/^2 



(45) 



M'2 



In the case of the zero recoil vicinity (r ^ 0) 



M2m'i ' 



(46) 



while the "common" factor, which may be considered as the 'universal' Isgur-Wise 
function, is written in the frame of the initial meson at rest ( p = ): 



i 



luj'^ni^nil 



iiPi -Pi) + mim[)ujiuj[uj2 , 



(p'p'iQp'T)-MH)^ 



(27r 



(47) 



where (pp'iq'p'T) and 0p(|q|) correspond to the radius parts of the wave functions of 
the meson in the initial state and the one in the final state, respectively. Being 
in a covariant form and in the C.M.S. one, they may be obtained by solving the 
Schrodinger equation Eq.(22) and its specific one in its center mass system respec- 
tively, as long as the QCD inspired potential is transformed into the corresponding 
form in the equation. We should note i). Of the approximations, the instantaneous 
one on the whole matrix element, is essential in the approach, and we think it is 
relible to the purposes due to the facts that we restrict to apply this approach only 



17 



to the weak binding systems and the approximation only for the B-S equation, is 
proven to work well in many cases as long as a weak binding system is concerned f^^^. 
ii). The wave functions obtained by the present way are more reliable than those by 
other ways, because the adopted potential is proven to work well for heavy quark 
systems, although the 'universal' function C, being the overlap integration of the 
wave functions, does not sensitive to the specific radius wave functions very much, 
especially when the momentum recoil is not great it is contralled by the normal- 
ization condition of the wave functions, i.e. the overlap integration is approaching 
to the normalization when the recoil is approaching to zero and the reduced mass 
of the final state is approaching to that of the initial one in the meantime, iii). 
When carrying through the trace for 7 matrices in the Eqs.(39) and (42) so as to 
reach to the form factos, the contributions from (7 ■ q) in the terms fi and ^i, have 
been ignored savely in the considered accuracy, because they are small in the case of 
weak binding, and when carrying through the integration Eq.(47), of the integrand 
all terms proportional to q in odd power will not contribute at all ( only even power 
terms contribute ). 

The proposed approach to calculate the weak decay matrix elements, so as to 
the form factors is completed. It is expected that the approach is available as long 
as the mesons in initial and final states are of weak binding. It is interesting to 
compare ours with that of IGSW model'^^'. In their model, the authors calculated 
the form factors by using the Gaussian-type wave functions, of which the parameters 
are determined by the variational method. It is easy to see that in the case of 
weak binding and at zero recoil vicinity, the formalism of our approach is consistent 
with theirs except a tiny difference in the formulation for the overlap integration 
of the wave functions'^^"^^' . However still in the case of weak binding but with a 
large recoil, there are two remarkable deviations between the two approaches. One 
deviation comes from the difference in the spin structure of the wave functions. The 



other from the arguments in the wave function integrand. For instance, the function 
corresponding to the ^ in IGSW model is \ n^^ ^3 ^^^ readst^^': 

3 77T.2 tm t 

i,asw = (-^) " exp"2MrK^(/3^ + /3-) , ,48) 

where tm = (M — M')^; M and M' are the masses of the "mock meson" f^'^'; P and 
f3' are the variational parameters for the initial and the final states respectively; 
7712 is the mass of the 'spectator' and /t is a parameter introduced by hand. In 
the IGSW model, k is adjusted to be 0.7 by fitting the n electromagnetism form 
factor and the authors of Ref. [17] regarded it as relativistic corrections due to a large 
recoil. However in our approach it is different i.e. all factors come into the formula 
automatically. The interesting thing is that occasionally the numerical calculation 
shows that the function of C, obtained in our approach, is very close to that of the 
IGSW model with k = 0.7. To show the fact, we present the corresponding ^ 
functions of B^ -^ J/ip + X and Be —^ Bs + X in Fig. 2 and Fig.3 respectively. The 
dotted line represents the ^ function obtained by the Gaussian wave functions and 
with the original formalism of Ref. [17] and /t = 1. as well. The dashed line represents 
that obtained by the wave functions from solving the Schrodinger equation of the 
potential model'^^^ with the same formalism. The dotted-dashed line represents 
that of the IGSW model but with k = 0.7. Note that the parameters /Sb^ = 0.88 
and f3j/qi = 0.65 which were obtained by the authors of Ref. [22] are used here for 
the dotted, dashed and dotted-dashed three curves. The solid line represents the ^ 
function achieved by our approach with the wave functions obtained by the potential 
model f^^l. It is easy to see from the Fig.2 that the result of Eq.(47), is very close 
to that of IGSW model with k = 0.7. It means that our results involve reasonable 
effects automatically. In order to have more comparison we also show the C, function 
for Be —* Bg + X in Fig.3 and the meaning of each line is the same as in Fig.2, 
although we do not expect that our approach is so suitable as that for the formal 
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decays Be -^ J/ip + X, due to the s quark is not so heavy. 

It is also very interesting to see the behavior in the hmit when the quark mass 
is approaching to infinite because it will let us see the relation between the form 
factors obtained here and the universal Isgur-Wise function. At the limit of 

1712 << nil, m[ and mi,m[^'00, (49) 

the formulae of Eqs. (44) - (47) reproduce those of the Isgur-Wise formalismt^"^' for 

the form factors. In fact, in the limitation of Eq.(49) and from Eq.(33) and Eq.(43), 

we have 

Pi — > mi ■ V, 

p[ — > m[ ■ v', (50) 

e' -> e'. 

Hence the Eq.(42) can be rewritten as 

Ur) = ^{v ■ v')tT [(1 + f){a + /575/)(l + ^)] , (51) 

where ^(f ■ v') is just the universal Isgur-Wise function and has the following form: 



Civ ■ v') = -^^ f yf^r (g;'T)0(|q|). (52) 

Thus the Eq.(51) and Eq.(52) reproduce those of the Isgur- Wise's formalism in the 
infinite heavy limit'^^^^, so a link between the Isgur-Wise function and the non- 
relativistic wave function overlap integration has been established. A factor of 
[0.5(1 + f • f')]^2 has been derived by Bjorken^^^l by using the Cabibbo-Radicati 
sum rule'^^], however Here in our formalism it automatically appears in Eq.(52). 

In Ref.[9], we have applied the approach to calculate the weak decays of the Be 
systematically. It is no doubt that when the experimental study of B^ meson has 
fruitful results, our approach will receive serious tests. We would emphasize here 
that for the decay modes B^ -^ J /'^{rje) + X, which might be used as the most 
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typical channels to identify the Be meson later on in experiments, the results from 
our approach are the most reliable ones, as in the calculation of the form factors the 
adopted nonrelativistic wave functions obtained from the potential model are tested 
well in addition to our approach itself is suitable them. Although our predictions 
to the processes are very close to those of IGSW model with k = 0.7'^°^, there 
exists some deviations in numerical results for the processes B^ ^ Bg + X, i,e, our 
predictions to the laters are larger than those of Ref. [20] both for semileptonic decays 
and nonleptonic decays. The deviations may be understood as follows: although the 
^ function for the modes is smaller than that of the IGSW model as shown in Fig. 3, 
the form factors gain an enhancement from the spinor factor as shown in Eq.(39) 
and Eq.(42). 

We should point out that the approach has quite wide potential applications, as 
long as the binding is weak . For example, the approach may be applied to those 
decays between the baryons which contain two heavy quarks in addition to the B^ 
meson's by means of a heavy diquark picture, and we will describe it in Ref. [26]. 

As stated above, when taking the instantaneous approximation to the whole 
matrix element by doing the integration of go on Eq.(23), an extra approximation 
has been adopted i.e. to take the expansion Eq.(26), thus the branch points are 
alternated. The extra approximation set down one more restraint 

^ « 0(1) 
nil 

i.e. the momentum recoil cannot be too great. If treating the branch points in the 

integrand by an accurate method i.e. making cuts, in principle, we may integrate out 

go exactly. However, the potential models only offer us the numerical solutions owing 

to the complicated concrete potential, thus we have to deal with an integration not 

analytically throughout, thus we cannot guarantee results so accurate by numerically 

integrating a counter with cuts, finally as a result the accurateness is lost too. We 

would conclude that the extra approximation may be weaken, but not be avoid at 
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all'^^l The experimental information on various decays is desired very much now, 
as we are at the position, that the approach remain to be tested and improved 
guiding by the information. In the near future the experiments especially on the 
Be decays are very important. Having the information, we will learn much on the 
decay mechanism as will as on the approach. 
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Figure Captions 

1. Feynman diagram corresponding to the weak current matrix element sand- 
wiched by the Be meson state as the initial state and a single particle state of 
the concerned final state. 

2. The ^ function for Be -^ J/ip + X. The dashed line: the ^ function of IGSW 
model with k = l.\ the dotted-dashed line: the ^ function obtained by the 
formalism of Ref.[18] with the wave function solved by potential l'^^^; the 
doted line: IGSW model with k = 0.7; the sohd line: the ^ function obtained 
by our approach. 

3. The ^ function for Bf. ^ B^ + X. The meaning of each type line is the same 
as that in Fig. 3. 
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